Chapitre 1 • Outils mathématiques de base 


On utilise souvent une définition alternative de ces intégrales généralisées dite au sens de la partie 
principale (ou valeur principale) de Cauchy 

/ oo /• R 

f(x)dx = lim / f(x)dx 

-oo 


PP [ 

J a 

Si l’intégrale converge au sens standard elle converge aussi en partie principale et les deux défi¬ 
nitions donnent la même valeur de l’intégrale. Une intégrale peut converger en PP sans converger 
au sens standard. 


f(x)dx = 



L 


f(x)dx+ / f(x)dx> si / non bornée en x = c G ]«,&[ 


b) Séries 

Séries numériques 

On appelle somme partielle de rang N d’une série numérique de terme général u„ la somme : 

N 

Stf ^ ^ Un 
n= 0 

La série converge et a pour somme S si la suite Sn a une limite bornée : S = lim .S’,v et on note : 

N —>oo 
oo 

S = ^ ^ U n 

n =0 

Une condition nécessaire de convergence est que lim u n = 0 ; la convergence de la série de terme 

n —>oo 

général |w„| entraîne celle de la série de terme général u n qui est alors qualifiée d’ absolument 
convergente. 

Exemples 

OO 1 

. soit £"(x) = — (fonction Ç de Riemann); cette fonction est définie (/'e la série converge) 

' n x 

n =1 

pour x > 1 ; 

oo 

2, la série géométrique de raison q '■'^q n est convergente pour jg| < 1 et a pour somme 

n=0 

i 

T - g' 

• Critères de convergence pour une série à termes positifs (u n ^ 0) : 

o critère de d’Alembert : soit L = lim ——. Si L < 1 la série converge; si L > 1 la série 

n —>oo u n 

diverge ; 
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1.1. Rappels d’analyse 


o critère de Cauchy : soit L = lim (n,,) 1 /". Si L < 1 la série converge; si L > 1 la série 

Tl —»00 

diverge ; 

o critère de comparaison : si lim — = L avec 0 < L < oo alors les deux séries de terme 

n kso v n 

général u n et v n sont de même nature. 

Pour une série alternée u n = (—1 ) n a n avec a n >0 on a le critère de Leibniz : si à partir d’un 
rang N la suite a n est monotone décroissante et lim a„ = 0 (condition nécessaire) alors la série 

n — kx ) 

converge. 


Séries de fonctions 

Le terme général est une fonction : x —> u„(x). De ce fait la somme de la série est aussi une fonc¬ 
tion x —+ S(x). À x fixé on est ramené à une série numérique. La convergence de la série est alors 
assujettie à la valeur de la variable (convergence simple). On parle de domaine de convergence. 

S’il existe une série numérique convergente à termes positifs a n telle qu’à partir d’un rang N on 
&\u n (x)\ f a n Vx g \a. b\ alors la série de terme général u „ (x) est uniformément (et absolument) 
convergente dans [a , b] et la somme de la série est une fonction continue dans cet intervalle (critère 
de Weierstrass). 


• Cas particulier : les séries entières u n (x) = a n x" 


Le domaine de convergence est l’intervalle ] — R, R[ où R 


lim 

n — kx ) 


a n 


O-n+X 


est appelé rayon de 


convergence de la série (dans C le domaine de convergence est un disque de rayon R). Dans 
le domaine de convergence la somme d’une série entière est une fonction continue. La série est 
dérivable et intégrable terme à terme. 


1 -x 


Exemples de séries entières 

oo 

E*" 

n=0 
00 v n 

Y — 

^ n! 

n =0 
00 v n 

E- 

^ n 


= e 


= -ln (1 - x) 


|x| < 1 

|x| < 00 

|X| < 1 


E^ 

r =0 

oo 

E^ 

n =0 


(2 n)\ 

x 2n+1 

( 2 n + 1 )! 


= cosx |x| < oo 


= sinx |x| < oo 


c) Équations différentielles 

Il existe des méthodes générales de résolution dans un petit nombre de cas rappelés ci-dessous : 


© 


3 











Chapitre 1 • Outils mathématiques de base 


• Équations du premier ordre F(x, y , y') — 0 

• les équations aux variables séparables : l’équation peut être mise sous la forme f(y)dy = g{x)dx 
et intégrée directement ; 

• les équations linéaires : elles sont de la forme 

ao(x)y' +a\(x)y = f(x ) 


la solution générale est de la forme : y{x) = a y \ (x ) + yo(x) où y i (x ) = exp 


ai(x) 

CI()(X) 


dx 


est une solution de l’équation sans second membre, a une constante arbitraire et yo(x) une 
solution particulière de l’équation complète, qui peut s’obtenir par variation de la constante : 

" /U 0 


yo(x) = yi(x) 


-dx'. 


ao(x')yi(x') 

• Équations linéaires du second ordre à coefficients constants 

a 0 y" +a\y' + a 2 y = f(x) 


L’équation sans second membre se résout en posant y(x) = e rx où r G M ou C est déterminé par 
identification. La solution particulière de l’équation complète s’obtient par identification selon la 
forme du second membre ou par la méthode de variation de la constante. 


1.2 Les fonctions utilisées en physique 

• Fonctions continues ou continues par morceaux, à dérivées continues ou 
continues par morceaux 

Une fonction est de classe C k sur un intervalle I (éventuellement M) si elle est continue sur I ainsi 
que ses k premières dérivées. Une fonction C°° est indéfiniment dérivable. 

Une fonction continue par morceaux a des discontinuités de première espèce (de saut fini) et est 
continue sur les intervalles délimités par les discontinuités. 

Exemples 

1, la fonction valeur absolue t —» \t\ est C 1 par morceaux car elle est continue (C°) et sa dérivée 
première est continue par morceaux; 

2, la fonction (ou échelon) de Heaviside 

u : t —> u(t) = 0pourf<0 
= 1 pour t> 0 
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est continue par morceaux (ou C° par morceaux); 

3 les fonctions causales sont nulles pour t < 0. Elles satisfont l’identité f = uf où u est la 
fonction de Heaviside. 
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1.2. Les fonctions utilisées en physique 


• Fonctions périodiques 

3 T > 0 : f{t + T) = fit) Vf. Le plus petit réel positif T qui satisfait cette identité est la 
période de la fonction qui est alors qualifiée de T -périodique. Exemple : la fonction harmonique 

t —> exp ^2/77—j . La fréquence est définie par v = — et la pulsation par co = 2ttv = 

• Fonctions caractérisées par leurs propriétés d’intégration 

• Les fonctions sommables (ou intégrables) sont celles dont l’intégrale / |/(f)| dt est bornée. 

J R 

Elles forment un espace vectoriel noté C 1 . 

• Les fonctions de carré sommable sont celles dont l’intégrale / |/(f)| 2 dt est bornée. Elles 

forment un espace vectoriel noté C?. Ces fonctions modélisent en particulier les signaux d’éner¬ 
gie finie ou encore les fonctions d’ondes de la mécanique quantique. 

• On peut définir plus généralement l’espace C k pour k quelconque dont les deux espaces C 1 et 
C 2 sont des cas particuliers très utilisés en physique et théorie du signal. 

• Les fonctions localement sommables sont celles dont l’intégrale J |/(f)| dt est bornée sur tout 

intervalle I fermé borné. Elles forment un espace vectoriel noté C) nc . On utilise parfois C 2 oc 
l’ensemble des fonctions de carré localement sommable. 

Il est clair que les fonctions de £} sont localement sommables (et celles de C 1 de carré locale¬ 
ment sommable). 

• Fonctions tempérées 

Ce sont les fonctions localement sommables et à croissance lente, c’est-à-dire celles qui croissent 
au plus comme une puissance (finie) de t. Autrement dit on peut trouver un entier N non négatif 

, / . | ,. fit) 

borne tel que 


lim 

Ifl—>oo t 1 


jN =0. Notons qu’une fonction décroissante à l’infini est forcément à 

croissance lente. 

La plupart des fonctions utilisées en physique sont des fonctions tempérées. 


• Égalité presque partout 

On dit que deux fonctions / et g sont égales presque partout et on note / = g si /(f) = g(t) 
Vf sauf sur un ensemble dénombrable de valeurs de t (ensemble de mesure nulle). Par exemple la 
dérivée de la fonction de Heaviside est presque partout nulle. Les fonctions égales presque partout 
ont comme propriété : 

f = Ji f(t)dt = J git)dt V/ 


Cette propriété entraîne que les éléments de C 1 (de C 2 ) sont en fait les classes d’équivalence des 
fonctions égales presque partout et (de carré) sommables. 
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• Modélisation des signaux impulsionnels 

L’impulsion <5 de Dirac (définition « phénoménologique » ; pour une théorie rigoureuse voir le 


chapitre 5). 


V/ 


r 

J — ! 


ô(t — a ) f(t) dt = f(a) 

«OO 

V/ / = (—1)" / (n) (0) 

J —co 

g(t)8(t - t 0 ) = g(t 0 )8(t - t 0 ) 

8(at ) = ,—- S(t ) en particulier 8(—t) = 8(t) 

\ a \ 

La notation Vf doit rester compatible avec l’existence des valeurs de f qui y apparaissent 
(/(O), f(a), /'(O),...). 

Si / est discontinue en t 0 avec pour saut cr on peut écrire fit) = f c (t) + cru{t — t 0 ) où f c est 
continue, alors : 

au sens des fonctions —> f'(t ) = f'(t) presque partout 
au sens des distributions —> f'(t) = f'(t) + cr8(t — to) 

En particulier, au sens des distributions if (f) = <5(f)- 

• Formule de Leibniz sur la dérivée n ieme d’un produit de deux fonctions 

ri 

{u (x) V (x)] (n) = J2 Cn [« (^)] ( " ) [W (x)\ (n ~ k) 


k =0 


où C„ est le coefficient du binôme : C„ = 


n ! 


kl (n — k)\ 


Définitions 


Fonction de Heaviside : u(t) 


Fonction « signe » : sign(f) 


0 pour t < 0 
1 pour t > 1 
— 1 pour t < 0 

1 pour t > 0 

2 u(t) — 1 


Fonction « porte » : I \r(t) = 1 pour t G 

= 0 sinon 


T T 

2 ’ 2 


Fonction « üiangle » : A T (t) 


Fonction sinus cardinal : sinc{t) 


(l - y) « V - l'I) 

sin(7r?) 

ut 
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1.3. Les Séries de Fourier 


1.3 Les Séries de Fourier 


/ est une fonction T -périodique sommable sur tout intervalle A fermé de longueur T ; on note 
27 T 

(o = — la pulsation fondamentale. On appelle coefficient de Fourier de / : 


Cn = ^ [ f{t) e ~ ina>t dt ne Z 
1 J A 


On définit la série de Fourier de f : 


S F {fU) = Y, c ne inwt 

nÇzJj 


o Si en f 0 e A, / et f'o nt une valeur à droite et une valeur à gauche alors (théorème de 
Dirichlet) 

S F (f-,t 0 )= \ [/(t 0 + ) + /('<,")] 

o Si / est continue et /'continue par morceaux alors S F (f',t) = fit) uniformément Vf et la 
série de Fourier de f 1 s’obtient en dérivant celle de / terme à terme. 

o Si / est C k ~ { et f {k) est continue par morceaux alors le „| — ——r pour n grand. 

n k+i 

• Série de Fourier de fonctions réelles : / réelle =>■ c_„ = cf, (complexe conjugué de c„) ; alors 


1 OO 

s F (fu) = +y, a n cosncot + b n sinnwt 


n =1 

ci n — Cn "F C—n b n — i ( c u C— n ) 'te c n — — (o n ib n ) C— n — — ia n + ib n ) 

2 f 2 f 2 f 

cio = — f(t)dt ci,, = — / fit) cos nojt dt b n = — / fit) sin ncot dt 

T Ja ^ J A T ,/ A 


On rappelle la relation de Parseval (pour f e C 2 (A)) 


f f\m\ 2 * = 5>i 2 = (t) 2+ ÎE(“- + ^) 

nez. n =1 


© 
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1.4 Les fonctions définies par des intégrales 

Quelques fonctions usuelles 


Fonctions eulériennes 


La fonction Gamma TQ) 


La fonction Bêta B(x,y ) 


p+OO 

r(.x) = / e ’t x ~' dt x € 

J o 


p+OO 

r <n) (x) = / e~‘(\nt) n t x ~ l dt x G 

J o 


T(;c + 1) = x T(x) Vt £ 8 \ {Z } 


r(l) = 1 , T(-) = yfïr 

r'(l) = -y r'(i) = —(y + ln4)v / 7r 
y = 0,577... constante d’Euler 


1 (2 nV 

r (« + !) = «! r ^ 7 + 2 } = V77 ' 


Pour Q, y) G R* +2 

B(x,y) = [ t x ~ l (1 - t) y ~ l dt 

J o 

2L 

= 2 2 (cos df x (sin d) 2> '~'d 0 

J o 

/•°° C” 1 

= / -dt 

y o (i+o* +y 


B(x,y)= V*,y,x + y G R\ {Z } 


77 


B(x, 1 — x) = r(jt)r(l — .r) = —- x ^ TL 


sin TTX 


Toutes les fonctions ci-dessous sont définies sur R* + ; certaines peuvent être prolongées sur 
ou même sur R. 


2 r ,2 

Fonction Erreur Erf(x) = —= / e dt 

J o 

2 f°° ,2 

Erreur complementaire Erfc(x) = —= / e dt 

Jx 

Sinus Intégral Si(x) = / dt 

J 0 f 

r + °° e~‘ 

Exp elle Intégrale E n E n (x) = J — dt n G N* 

f +co cos t 

Cosinus Intégral Ci(x) = — - dt 

J x 1 

r e 1 

Exp elle Intégrale E, E; (x) = PP — dt 

J — OO 1 
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Énoncés des exercices 


Produit de convolution 


f *g :t ->(f * g)(t ) 
Si / et g sont causaux : (/ * g)(t) 





On utilise couramment la notation (abusive) f(t) * g(t). 


t') g(t') dt' 

- t') g(t') dt' 


• Propriétés 

• symétrie '■ f * g = g * f 

• si / G £* et g G C) alors / * g G C) 

• si f £ C 1 et g' est continue et bornée à l’infini alors (f * g)' = f * g' 


ENONCES DES EXERCICES 


© 


1.1 Intégrales généralisées 

D Étudier la convergence des intégrales suivantes : 


(a): 

(c) 


1 


0 yjx{\ - X) 

X + 1 


dx 


(b): 


1 


l-i (x + l)J/x 


dx 


r 

(d): / c~ x cosx dx 

J o 


/*oo 

•'° \J'x (x 3 + l) 

x a ~ l 

Etudier en fonction de a la convergence de l’intégrale / - 

J o * + 1 


dx 


L’intégrale 


1 


t 3 + 1 


dx est-elle convergente au sens standard ? au sens de la partie princi¬ 


pale ? Si oui la calculer. 

1.2 Séries 

D Déterminer la nature (convergente ou divergente) des séries : 

OO - OO 1 OO , OO 

<a) : £^j= (d>:£ 

“ Vn 2 + 1 “ 3" + n n 2 ^ 

n =0 n=0 n =1 n=0 


(-i r 


Pour quelles valeurs de a la série de terme général (1 — a)" est-elle convergente ? Calculer sa 
somme. 

OO 


El Calculer la somme de la série V'' —. 

' n\ 


n =0 
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Q Soit deux séries numériques absolument convergentes de terme général respectivement a n et 
b n . Montrer que la série de terme général a„ cos nx + b n sin nx est uniformément convergente 

Vx G M. 

1.3 Équations différentielles 

D Résoudre les équations différentielles du premier ordre : 

(a) y' + y 2 sinx = 0 avec y(0) = 1 (b) y' = a 2 — y 1 (c) (1 + x) y 1 + - = 


x x z 


0 Résoudre les équations différentielles du second ordre : 

(a) y" + 2y' + 2y = 2x — sinx (b) y" + y = ln |cosx| 

1.4 Impulsion de Dirac 

D Établir les propriétés suivantes de l’impulsion de Dirac : 

(a) / f(t)S(t — t 0 )dt = f(t 0 ) si t 0 G]n, b[ et 0 sinon. 


(b) 8{at + b) = - — -S(t H—) 
n a 


Calculer : 


(a) 


(2 1 + 3 )8(t + 1 )dt (b) 


8{t + 1 )dt 


(c) J (t 3 + 3 t 2 )8'(t - 1 )dt (d) 


/»oo 

f te~ ,7 

J o 

/ OO 

t cos (t) 8 { "\t)dt 

-OO 


0 Calculer t"8 ik \t). En déduire n solutions indépendantes au sens des distributions des équa¬ 
tions : 

t"T = 0 
t"T = S 


£) Potentiel coulombien à une dimension : 

a) Calculer « au sens des distributions » la dérivée seconde de x i—► |x|. 

b) On place une particule ponctuelle de charge q à l’origine (x = 0) dans un espace à une 
dimension (une ligne). Le potentiel créé par la charge à l’abscisse x est solution de l’équation 
de Poisson à une dimension 

V"(x) = ~ — S(x) 

Trouver la solution de cette équation qui soit une fonction paire. 
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Énoncés des exercices 


c) On place sur la ligne une autre particule, de charge —q. À l’instant t = 0 elle est située en 
xq ^ 0 et possède une vitesse vq. Montrer que quelle que soit sa vitesse initale vq elle ne peut 
pas échapper à l’attraction de la particule localisée en x = 0. 

On rappelle que la force que ressent en x la particule de charge —q plongée dans le potentiel 
V(jc) est : F(x) = qV'(x). 

1.5 Séries de Fourier; propriétés, dérivation 

D Développer en série de Fourier la fonction T -périodique suivante : 

f{t) = at -T/2^t^0 

= bt 0^t<T/2 

où \a\ ^ \b\. Discuter les cas a = ±b. 

0 Développer en série de Fourier la fonction T -périodique suivante : 

g{t) = a —T/2 ^ t < 0 

= b 0^t<T/2 


0 Comment peut-on comparer ce résultat à celui de la question précédente ? 


1.6 Série de Fourier et somme de séries 


Développer en série de Fourier la fonction paire et 277-périodique définie pour t G [0, ir] par 
f{t) = e~ al . 

En déduire la valeur de la somme des séries numériques : 




n= 0 


(- 1 )” 
n 2 + a 2 




n =0 


n 2 + a 2 


1.7 Série de Fourier 

Développer en série de Fourier la fonction 277-périodique définie par : f{t) = sinfroQ?) pour 
—tt ^ t < tt avec cûç) G ]0,1[. Que se passe-t-il quand (o 0 —> 1 ? 

1.8 Série de Fourier et équations différentielles 

Trouver une solution particulière périodique de l’équation y"(t) — 4 y(t) = f(t) où / est 
27r-périodique définie par f{t) = |f| pour t G [—77, 77 ]. 

1.9 Séries de Fourier sinus et séries de Fourier cosinus 

Soit / une fonction satisfaisant les critères de Dirichlet sur ]0, L\. Construire une série de Fourier 
de périodes 2 L ne comportant que des termes en sinus (cosinus) et qui coïncide avec / sur ]0, L\. 
Applications : 

• Développer en série de Fourier-cosinus sur [ 0 , 77] la fonction f{t) = sin t. 

• Construire la série de Fourier-sinus puis la série de Fourier-cosinus qui coïncident avec la fonc- 
© tion t —> coshï sur ]0,1[. 
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1.10 Fonctions eulériennes 

f) Calculer r(-), B(2, — -), T(« + -) et en déduire T (—n + -) (n G N*). 

/»00 

0: Calculer J n = / t"e~ at Int dt (n G N, a G M* + ) en termes de T et de T'. 

J o 

Donner la valeur de / () et de J \. 


1.11 Calcul d’intégrales 

Calculer les intégrales suivantes : 


h = [ t a e~^dt i/3 > 0, a > —1) ^ I 2 = f 

J o J o 


-dt (a, b >0) 


D h 


f 1 

/o \A (l — t) 

r*7r/2 


dt 


y/t + b 2 

/•OO 

E> / 4 ( jc ) = / e~ ia,2+2b,+c) dt (a > 0) 

J X 


[‘77/ £ /*00 

h = (tan 9) a d9 a G ] — 1,1[ 0 4 = / cos(pt)Ci(t)dt 

J o J o 


h — 


j 

f c 

D i,-J 


t + b 
\/lnf 


dt (a,è > 0) 
dt 


D /s - 

ES Ao 




t + b 


dt {a > 0 ; -b f. [a, /3]) 


smaf 
t + b 


dt (b > 0) 


/n = 


,°o J 

./() 1 + f 2 " 


dt {n g N*) 


/ 12 = / (9 - t 2 )ï dt 


1.12 Produit de convolution 

D Calculer les produits de convolution / * g (les fonctions 1I ( , et A a ci-dessous sont définies 
dans le rappel de cours) : 


(a) fit) = Il 2 {t) 

(b) fit) = n 2 (r) 

(c) fit) = sint 

1 


g(0 = /(0 

g(t) = 2A 2 (t) = (2 - |î|) n(2 - |f|) 

git) = aA a it) = ia — |f|) nia — |f|) a > 0 

1 


(d) fait) = TT--t a ‘c ’uit) a > 0 fbit) = -J— t b ] e r n(t) b > 0 

lia) tib) 


Calculer / * 
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Énoncés des problèmes 


§) On désigne par f a la fonction définie pour a > 0 par f a (t) = — i^ e 

ClyJlTT 

a) Montrer que / '* — 1 ^ 


— t 1 /2a 1 


r 

J —( 


+ P'dt = A / — exp — 
a V 4a 


b) Calculer f a * f h . 

c) Calculer f* 3 = f a * f a * f a puis f*" (n entier positif). 

ÉNONCÉS DES PROBLÈMES 


Problème 1.1 Équation de la chaleur 

Trouver à l’aide des séries de Fourier la solution du problème suivant : 


définie pour x G [0, l] et t > 0 
= a-^-r — fiu a > 0 ; P > 0 


u(x,t ) 
du 

dt dx 2 

u(0,t) = u(l,t) = 0 t^O 

u(x, 0) = x(I — x) xE[0,/] 


Problème 1.2 Séries de Fourier, polynômes de Bernouilli et fonction f 
de Riemann 

Les polynômes de Bernouilli sont définis par : 


B\(x) = x — 


1 


B'(x) = 

B n ( 1) = B n ( 0) 


« > 2 


On veut déterminer la série de Fourier S n de période 1, qui coïncide avec B n pour x E ]0,1[. 

Q Déterminer 

0 Établir que le coefficient a 0 de S n est nul Vn. 

0 Quelle relation existe-t-il entre S n et S n -\ ? En déduire l’expression de S n ; on distinguera le 
cas n pair du cas n impair. 

Q Justifier que S n ( 0) = B n { 0) pour n ^ 2. Calculer B^{x) et Bpx) et en déduire /(2) et /(4) où 
la fonction / de Riemann est définie par : 


© 


fw = E 

n =1 


î 

n x 
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Problème 1.3 Phénomène de Gibbs 

D Soit / une fonction 277--périodique admettant un développement en série de Fourier et S N (t) 
la somme partielle d’ordre N de la série de Fourier 

N 

Sn(.ï) = flo /2 + ^2 a n cos {lit) + b n sinfnf) 

n= 0 


a) Montrer que 


b) En déduire que 


a n cos fit) + b n sin(/îf) = 


77 


f(t')cos(n(t' — t))dt' 


S N {t) = 


277 


sin \(N + \){t' — t) 1 
fit ')— 1 r , 2 —^-Udf 


sin 


[\{f - 0 ] 


Soit la fonction /, impaire, 277-périodique, définie par fit) = 1 sur [0, tt[. 

a) Déterminer sa série de Fourier. 

b) En utilisant le résultat de la question Q b) montrer que pour t G ] — 77, 77 [ 


Sivit) — — 


rt sin [jN + \)u] _ J_ 

_ f sin(i M ) 277, 


sin [(A^ + ^)n] 
sin (\u) 


du 


c) On regarde le voisinage de t = 0 et on se place pour N grand. Montrer que dans cette 
situation la seconde intégrale est négligeable devant la première et que 


S N it) ~ 


1 


sin v 


77 m J 0 sin (£) 


dv 


où on a posé m = N + -,v = um. 


77 


d) En déduire que le premier maximum de SV(t) apparaît pour t ~ — ~ 

m 

valeur approchée du premier maximum de S^(t) est donnée par 


77 

N 


. Montrer qu’une 


oinax ^ 

— 


2 

77 


sin f 


df 


Quelle remarque peut-on faire au vu de ce résultat ? 

Problème 1.4 Formule de Stirling pour la fonction r 

D Montrer que 


( lnr(f)r/f = ln \/277 

J 0 
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Du mal à démarrer ? 


On définit pour x ^ 0 


r x+\ 


F (x) = 


ln T (t)dt 


a) Calculer F'{x). 

b) Al’ aide du résultat de la question D en déduire l’expression de F (x). 

0 En appliquant la règle des trapèzes (voir formulaire) à l’intégrale qui définit la fonction F 
montrer que (formule de Stirling) 


i» 1=> I/x + 1) = \Z2tt x x+1 e x (l + o (l /x)) 


Formulaire 


lnsin(7 Td)dO = — ln2 


0 < —-lnl/x) ^ -- pourx > 1 

dx- x — 1 


f 

J a 


f{x)dx = - {b-a) [/(a) + f(b )] - — (b - a) 3 f"(9) 


avec 8 £]a,b[ 


■v Æ 


DU MAL A DEMARRER? 


© 


1.1 O et 0 Utiliser les critères de convergence. 

E) Décomposer —-en éléments simples. Dans cette même question pour intégrer les élé- 

x i + 1 

ments simples prendre la partie principale à l’infini (justifier). 

1.2 D et 0 Utiliser les critères de convergence. 

Ë) Utiliser la table de séries entières. 

1.3 0 Seconde équation : utiliser la méthode de variation des constantes. 

1.4 D a) Définir une fonction g(t) = f(t ) pour t E ]a,b\ et g(t) = 0 sinon. 

0 d) et D : Utiliser la formule de Leibniz. Q a) : La dérivée d’une fonction discontinue 
est différente selon qu’on la calcule au sens des fonctions ou au sens des distributions. 

1.5 Vérifier les symétries éventuelles de la fonction selon les valeurs de a et b. 

0 On remarque que g = f'. 
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1.6 La fonction est continue donc elle est égale à sa série de Fourier en particulier en t = 0 et 

t = TT. 


1.7 Prendre soigneusement la limite coq —> 1. 

1.8 L’équation différentielle ne fait intervenir que des dérivées paires et le second membre est une 
fonction paire ; on cherche donc la solution sous la forme d’une série de Fourier paire. 


1.9 Prolonger la fonction / sur ] — L, 0[ de façon à obtenir une fonction paire (impaire) et prolon¬ 
ger la fonction ainsi définie ]—L, L[ sur M par périodicité de période 2 L. 

1.10 Attention (voir table dans les rappels de cours) : on définit r (x) pour X < 0 non entier en 

r (x + /t) 

prolongeant les valeurs pour x > 0 à l’aide de l’identité T ( x ) =-pour 


x G \—n. — n + 1[ avec n G N*. 

1.11 I 2 : poser u = Vt + b 2 ; 7g : intégrer par parties et utiliser / 

J o 


x (x + 1)■ 
sin(at) 


■ {x + n — 1) 


77 


dî = — sign(a); / 8 : 

distinguer le cas a < (3 < —b du cas — b < a < (3 \ I 9 : poser t = e u ; I\\ : poser u = t 2n ; 
1 [2 : poser u = t 2 / 9. 


1.1 2 D a) et b) On peut utiliser un raisonnement graphique. 


Problème 1.1 

Les conditions aux limites suggèrent de chercher des solutions de la forme 

E X 

b„(t) sin/î7r— ; 

n^l 

pour implémenter la condition initiale utiliser les séries de Fourier sinus (exercice 1.9). 

Problème 1.2 

Remarquer que la série de Fourier é>i est une fonction impaire. Dans la suite, attention à ne 
pas faire de confusion entre l’indice n qui numérote les polynômes de Bernouilli et l’indice k qui 
repère les coefficients de Fourier. 

0 La relation de récurrence entre les cik(n ) et les bk(n + 1) sera transformée en relations par pas 
de deux dans une même série de coefficients. 


Problème 1.3 

0 b) Dans l’expression de Sn(1 ) on fera un changement de variable u = t — t' dans l’intégrale 
entre —77 et 0 et u = t' — t dans l’autre intégrale. 

0 c) On établira que pour t ~ 0 les deux intégrales de la question O sont de la forme 



2eg{a) + o{e). 


Problème 1.4 

D On fera le changement de variable t' = 1 — t. Dans tout le problème, utiliser les propriétés de 
la fonction T. 
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Corrigés des exercices 


Corrigés des exercices 


1.1 Convergence des intégrales : 

a) La fonction — est non bornée en 0 et 1 qui sont les deux bornes de l’intégrale. 


y/x(l - x ) 
Pour x ~ 0 : 


1 


1 


—— = —— ; le critère de comparaison avec l’intégrale de 

.. / V V 1 / ^ 


y/x{\ - x) y/x ^ X '/2 
Riemann est applicable avec a = — < 1 donc l’intégrale converge à la borne inférieure. 

Pour x ~ 1 : — ~ - = -— ; là encore le critère de comparaison avec 

y/x(\-x) s/T^ (l-x) 1/2 

l’intégrale de Riemann est applicable avec a = - < 1 donc l’intégrale converge aussi à la 
borne supérieure ; elle est donc convergente. 

b) La fonction -- est non bornée en — 1 et 0. Pour x ~ 0 : 


(x + l)^/x 


(x + l)ÿx 


: l/3 ’ 


1 


le critère de comparaison avec l’intégrale de Riemann est applicable avec a = — < 1 donc 

l’intégrale converge en 0. Pour x ~ — 1 : -— ~-; le critère de comparaison 

(x + \)\/x —(x + 1) 

avec l’intégrale de Riemann est applicable avec a = 1 donc l’intégrale diverge à la borne 
inférieure ; elle est donc divergente. 

c) Pour x rv Q \ - — -— ; on a donc de nouveau a = - < 1 et l’intégrale converge 


(* 3 + 0 

V X + 1 1 

en 0. À l’infini — - ~ —. On a ici a = 1 donc l’intégrale diverge à l’infini ; elle est 

\Ao+i) 1 

donc divergente. 

d) On a je x cos v| ^ e x qui est intégrable entre 0 et l’infini. Donc l’intégrable est conver¬ 
gente. 

°o v a -1 


- f x a 1 

Etudions la convergence de l’intégrale / - dx : en x ~ 0 : 

J o x + l 


r GL— 1 


x + 1 


v a ~ I 

~ X = 


A-a ' 


Le critère de comparaison avec l’intégrale de Riemann conduit à 1 — a < 1 
x "- 1 1 

À l’infini 


a > 0. 

1 _ a . Le critère de comparaison avec l’intégrale de Riemann conduit à 


© 


x + l x 

2 — a>l=>a<l. Finalement l’intégrale converge pour a G ]0,1[. 


17 


























Chapitre 1 • Outils mathématiques de base 


§) À l’infini 


1 


1 


1 


—- qui décroît plus vite que —. L’intégrale converge donc à l’infini. La 


x 3 + 1 

fonction ——- est non bornée enr = — 1. Au voisinage de X = — 1 : 


x 3 + 1 


1 


1 


1 


1 


x 3 + 1 


avec a = 1. L’intégrale / , 

J-oo * 3 + 1 

Au sens de la partie principale : 


(x + 1) (x 2 — x + l) (x + 1) (x + 1)" 

1 


dx est donc divergente au sens standard. 


1 


PP 


x 3 + l 
1 


(x + 1) ( x 2 — X + l) 


1 


x — 2 


i .-dx = - 

x 3 + 1 3 


PP 


x + l 


-dx — PP 


x + 1 x 2 — x + 1 

f°° X-2 , 

-dx 


x 2 — X + 1 


L’intégrale originelle est convergente à l’infini mais dans la décomposition de la fraction ration¬ 
nelle on fait apparaître deux intégrales, chacune étant divergente au sens standard. Comme 
l’intégrale de départ est convergente (à l’infini) on peut donc aussi bien la prendre en partie 
principale et donc chacune des deux intégrales ci-dessus sera calculée en partie principale à 
l’infini ce qui permet de récupérer la partie finie de la somme des deux intégrales. En outre, la 
première sera calculée aussi en partie principale en — 1. 


PP 


- dx = lim lim 

X + 1 R—>oo e —>0 


p-l-e 


1 


r -R 


'-R 


X + 1 


-dx + 


1 


J- l+e*+ 1 

= lim lim { In \x + l||ll -e + ln \x + 1||^, 

R—>00 €— A e 


dx 

y 


R 1 

= lim lim {ln e — ln (R — 1) + ln(î? + 1) — ln e} = lim ln-= 0 

R^ooe-0 1 \ J \ J 1 R l 

Dans l’autre intégrale, le dénominateur ne s’annule pas. Faisons le changement de variable 
1 

u = x — - : 

2 


PP 


F 


x — 2 
x 2 — x + 1 


-dx = PP 


r “ « - ! 


,2- 3 

,.R 


du 


' — oo W + ^ 


lim 

R—>oo 


l-R U 2 + | 


du — — 


3 2 2 u 

-= arctan —= 

2 y/3 V3 


2 J-ru 2 + \ 
= —77\/3 


du 


où la première des deux intégrales dans la deuxième ligne est nulle par symétrie (intégrale 
d’une fonction impaire sur un intervalle symétrique par rapport à l’origine). En regroupant les 
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Corrigés des exercices 


résultats on obtient 


1 TT 

PP / —; ~dx = 

. 7-00 x 3 + l V3 


Ce qu’il faut retenir de cet exercice 

Ce dernier exercice met bien en évidence le rôle de la limite symétrique dans la définition 
d’une intégrale en partie principale puisque la compensation des infinis se fait grâce à cette 
symétrie. 


1.2 9 Déterminons la nature des séries : 

a) Le terme général u n = t ~ — pour n grand. La série est de même nature que celle 

Vn 2 + 1 n 

de terme général - c’est-à-dire divergente (série de Riemann avec a = 1). 
n 


b) Le terme général u n = 


2 " + 1 

3" + n 


2\ n 1 + 1 / 2 " 


pour n grand. La série est 


1 + n /3" 

f 2\ n 

de même nature que celle de terme général - ) qui est une série géométrique de raison 
|<1. La série est donc convergente. 

c) Asymptotiquement le logarithme est dominé par la fonction puissance ; il existe donc 

e G ]0, 1 [ et N(e) tels que n ^ N(e) =£• ln n sC n e . On a alors — y sC —r—. La comparaison 

n 2 n l ~ e 

avec la fonction de Riemann pour x = 2 — e > 1 permet de conclure à la convergence de la 
série. 

(-!)" , , 1 

d) La série de terme général u n = —-=- est une série alternée avec \u n \ = —j= qui est mono- 

\ n Jn 


tone décroissant et tend vers zéro. Le critère de Leibniz est satisfait et la sérié converge. 

La série de terme général (1 — a) n est une série géométrique. Elle converge pour 11 — a\ < 1 

c’est-à-dire pour a G ]0,2[. Sa somme est-= —. 

1 — (1 — a) a 


1 


n =0 


n =0 


x ___ ' 

E) O n P eut écrire > —- = > — 

A—' n ! ' n ! 


= e l = 2,718... 


x=l 


O Les séries \a n \ et \b„\ convergent. Il en est de même de la série \a n + \b n \. Par 

n n n 

ailleurs \a n cos nx + h„ sin nx \ ^ \a„ \ + \b„ \ Vjc. Donc d’après le théorème de Weirstrass 

n n 

la série trigonométrique converge alors uniformément pour tout x et sa somme est une fonction 
continue. 

Ce qu’il faut retenir de cet exercice 

Cet exercice fait apparaître l’utilité du critère de comparaison avec la série de Riemann ou 
© avec la série géométrique. 
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1.3 Q On résout les équations différentielles du premier ordre : 
a) C’est une équation séparable : 

/ 7 y i 

y + y sin jc = 0 = — sin jc -= cos jc + C =>■ y(x) = 


-1 


y 


y(0) = 


-î 
l + C 


= 1 =>• C = -2 


finalement 


y(x) = 


y 


î 


b) Là aussi on a une équation séparable : 


2 — cos x 


cos jc + C 


/ 2 2 
y =a - y O 


y 


a 2 — y 2 


1 

2 a 


y + y 


a — y a + y 
a + y 


= 1 


= 1 => — ln 

2 a 

„2 ax 


a + y 


a — y 


a — y 


= Ke zax => y(x) = a 


= x + C 

Ke 2ax - 1 
Ke 2ax + 1 


c) On a une équation linéaire. On cherche d’abord la solution générale de l’équation sans second 
membre 

(i+ x) y + - = 0 

x 

qui est séparable 


y_ 

y 


î 


î î 


JC (1 + x) X 1 + JC 

X + 1 

yt(jc) = K - 


Recherchons maintenant une solution particulière de l’équation complète par la méthode de 
variation de la constante : 


>’oM = K (je) 


jc + 1 


(1 +x)y'o+— = ~ 2 ^ K ' (x) 

X x z 

=> K' (x) = 


1 


(x + l) 2 


1 


V- 

1 1 


JC (x + 1)" 


X + 1 (je + 1) 


Ce qui donne finalement 


jc 1 

K (x) = ln-+-+ C 

x + 1 jc + 1 


JC+1 JC 1 JC + 1 

y{x) =-ln-- 1 - C- 

jc jc + 1 JC JC 
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Corrigés des exercices 


0 Résolution d’équations différentielles du second ordre, 
a) On résout l’équation sans second membre en cherchant des solutions en e' x 

r 2 + 2r + 2 = 0=>r = -l±i 

yo(x) = e~ x (A cos * + B sin x) 

Le second membre est de la forme f\(x) + f 2 (x) avec (\ (x) = 2x et f 2 (x) = — sin x. On va 
donc chercher y\ et y 2 solutions de 

y" + 2 y[ + 2y\ = 2x 

y 2 + 2y' 2 + 2y 2 = — sin .r 

et la solution recherchée de l’équation complète sera la somme vi + y 2 . 

Cherchons yi sous la forme d’un polynôme de degré 1 : viOv) = ax + b. En identifiant on 
obtient : 

2a + ax + b = 2x ^ a = 2 , b = —4 =>- yi(^:) = 2x — 4 
Cherchons y 2 sous la forme d’une combinaison de fonctions trigonométriques : 

y 2 (x) = a cos x + b sin x. 

En identifiant on obtient : 


( 1 a + 2b) cos x + (b — 2a) sin x = — sin x => a + 2b = 0 , b — 2a = — 1 


2 1 

On en déduit a = - et b = — - et 


y 2 (x) = -(2cos;r — sin.r) 


Finalement la solution générale est donnée par 


1 


y(x) = e x (Acosx + B sin a:) + 2x — 4 + - (2cos^: — sin^:) 


b) L’équation sans second membre a pour solution y\(x) = A cos.r + B sin*. Pour trouver la 
solution particulière de l’équation complète on va utiliser la variation des constantes : on 
cherche yo(x) sous la forme 


© 


y o(jc) = A (.r)cosx + B (^)sin^: 

yb(x) = — A (^:)sin.r + R (,r)cosx + A'(^)cos^: + R'(.r) sin.r 

On a deux fonctions inconnues A (x) et B (x), il faut donc deux équations pour les détermi¬ 
ner. L’une sera l’équation différentielle ; pour l’autre imposons 

A' ( x)cosx + B' (x) sinx = 0 (1.1) 


21 


Chapitre 1 • Outils mathématiques de base 


ce qui simplifie l’expression de y' 0 (x) 

y' 0 (x) = —A (jc) sinjc + B (jc)cosjc 

Il vient donc pour y^(x) 

y[l(x) = — A (jc)cosjc — B (jc) sinjc — A' (jc) sinjc + B' (jc)cosjc 
En injectant dans l’équation complète on obtient 

—A' (jc)sinjc + B' (jc)cosjc = ln | cos jc | 

On a donc deux équations (1.1) et (1.2) pour A' et B' ; on en tire 

A'(x) = — sinjc ln|cosjc| ; B 1 (x) = cosjc ln|cosjc| 

On en déduit 


/ 


— sinjc ln |cos jc| d.r = 

U =COS JC 


A(x) = 

= cosjc ln | cos jc | — cosjc 

Pour calculer la primitive qui définit B(x) on intègre par parties : 


ln\u\ du = u ln I u I — u 


B(x) = 

= sinx ln I cos jc I + 


cosx ln | cos jc | dx = sinx ln|cosx| + 

1 


sin 2 x 


-dx 


cos.r 


cos.r 


-dx — / cos xdx 


On calcule 


- dx = 

COS JC t = Unx /2 


f l+t 2 : 

J T~^ï 


2 dt 

7 ? 


î 


î 


= ln 


1 — t l+t 

1 + tanjc/2 


dt 


Finalement 


B{ x) = sinjc ln | cos jc | — sinjc + ln 


1 — tanx/2 
1 + tan jc / 2 


1 — tan.r/2 

On reconstitue la solution >’q(x) = A (x) cosx + B (x) sinjc et on obtient : 

1 + tan jc /2 


y(x) = A cos x + B sin jc + ln |cos jc | — 1 + sin jc ln 


1 — tanjc/2 


(1.2) 


Ce qu’il faut retenir de cet exercice 

La méthode de la variation des constantes est peu connue pour les équations du second ordre. 
Dans ce dernier calcul, elle est indispensable. On peut aussi trouver une solution particulière 
à l’aide des séries de Fourier, mais sous une forme non explicite. 
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Corrigés des exercices 


1.4 [) a) On introduit une fonction auxiliaire : g(t ) = f(t) pour t G ]fl, b\ et 0 sinon. Ensuite 
on utilise les définitions 


8 (t — to)f(t)dt = 


/ 


8(t - to)g(t)dt = g(t 0 ) 
f(t 0 ) pour to G ]ci,b[ et 0 sinon 


1 b 

b) Montrons que S(at + b) = — <5(7 H — ) 

\ci\ ci 

f°° f°° t' dt' 1 -b 

/ 8(at + b)f(t)dt = sign(ci) / 8{t' +b)f (-)— = — /(-) 

./_™ t'=at a a fl fl 


r °° 1 b 

î8(t + —)f(t)dt 


a 


a 


Le facteur sign(a) devant la seconde intégrale vient du fait que les bornes ±oo sont inversées 

dans le changement de variable si a est négatif. Par ailleurs = —_. 

ci \a\ 

On utilise les définitions et le résultat de la question précédente l(a). 


a) 

b) 


/ CO 

(2 7 + 3 )8(t + 1 )dt = (2 1 + 3)| /= _, = 1 

-CO 

l 


te 1 S(t + 1 )dt = 0 car —1 n’est pas dans l’intervalle d’intégration. 


c) J (t 3 + 3t 2 )8\t - \)dt = - (7 3 + 37 2 )'| ï=] = -9 

r°° c i" 

d) / t cos(7) 8 (n) {t)dt = (— 1)" -[7cos(7)] 

7-oo dt n 


= n (—l) [n / 2]+1 si n est impair, 0 sinon. 


E) Calcul det"8 {k) (t ) : 

/ °o jk k jj 

( 0 = (-!)* E C i'â7['"]*î^7 [/(')] 




n ! 


7=0 


(n ~ j)l 


t"~J 


[/(Ol 


7=0 
(k-j) 


t=0 


t =0 


© 


Seul le terme constant (/ = 77) contribue dans la somme prise en 7 = 0. Comme j ^ k, si 77 > k 
il n’y aura pas de terme en t° et donc la somme est nulle. Si n ^ k alors 

/»oo /»CO 

/ t" 8 (k \t)f(t)dt = (-1 )*C* n 77 ! [/(7)]f = - 0 n) = (—1)"C£ 77! / 8 {k ~ n \t) f (t)dt 
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On a donc 


t n 8 (k \t) = (-\) n Cïn\8 {k ~ n) sin^k 
= 0 si n > k 

On en déduit que 

T = 'j2a k 8 (k \t) 

k =0 

(où les a k sont des constantes arbitraires) est solution de l’équation t n T =0 et 


n— 1 

T = (-1 ) n 8 (n \t) + l a k 8 {k \t ) 
k =0 

(où les a k sont des constantes arbitraires) est solution de l’équation t" T = 8. 

O Potentiel coulombien à une dimension. 

a) Calculons la dérivée au sens des distributions de la fonction valeur absolue : / : x i—> \x | ; 
f'(x) = sign(x ) ; f"(x) = 28 puisqu’il y a une discontinuité de saut 2 en t = 0. 

d 2 V q 

b) On déduit de la question précédente que la solution de l’équation de Poisson — y =-S(x) 

dx l e 0 

est : 

V(x) = ——\x \ + Ax + B 
2 e 0 

Pour avoir une solution symétrique x <-> — x (fonction paire) il faut A = 0. De plus on prend 
B = 0 ce qui revient à fixer l’origine des potentiels. 

c) Cette particule est soumise, de la part de celle qui est placée à l’origine, à une force 

. q 2 

F = qV =-sign(.r). L’équation du mouvement est donc 

2 e 0 


q~ 

mx = — - —sign(x) 

2e 0 

q 2 i 

Supposons qu’à l’instant t = 0 on ait x 0 > 0 et r>o > 0 alors x(t) = — - t~ + v qî+xq. Quel 

2m 6o 

que soit Vq il y aura toujours un instant où x(t) reviendra vers zéro. La force changera alors 
de signe et la particule sera toujours attirée vers le centre et ne pourra jamais s’échapper à 
l’infini. Ce phénomène de « confinement » est spécifique de l’espace à une dimension et est 
dû au fait que la force est constante. 

Ce qu’il faut retenir de cet exercice 

Ces exercices permettent de se familiariser avec les règles de calcul mettant enjeu l’impulsion 
de Dirac. La théorie rigoureuse sous-jacente à ces règles sera développée au chapitre 5. 
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Corrigés des exercices 


1.5 POn détermine les coefficients de la série de Fourier de / d’après les définitions et à 
l’issue d’un calcul simple d’intégrales on obtient : 

77 (a -b) [î - (-i)"] , (a + b) (~\r 

a 0 = -—(a - b) ; a n =-~- ; b n = - 

2w 7 tco n- (o n 


/(0 = ~£7 (fl_Z7)+ Ë 


(fl - b) [1 - (—1)"] 


cos ncot — 


(a+b)(-\) n . 


sin ncot 


n =1 


77 " CO 


O) 


Si |ût| y \b\ la fonction est discontinue et n’a pas de propriété de parité : (a n et/ou b n ) ~ —. 

1 11 

Si a = b la fonction est impaire et discontinue : a n = 0 et b n ~ —. Si a = —b la fonction est 

1 U 
paire et continue : b n = 0 et a n ~ —. 

n- 

0 Les coefficients de la série de Fourier de g sont 


a' 0 = (a + b) ; a' n = 0 


b' 


(fl-6) [!-(-!)”] 

77 n 


g(t ) 


+ b) 


(i a — b) 

77 


E 


[i 


(-in 

n 


sin ncot 


§) On a g(t) = f’(t). Si / continue, alors a' n = ncob n et b' n = —ncoa n par dérivation de la série 
terme à terme. Cette dernière relation pour b' n est ici bien vérifiée ; cela ne marche pas pour a' n 
car / est discontinue (si a / —b). En revanche 


, _ 2 

fl» j* 

_ 2 
"fl 


2 f ^ çT /2 f 

f\t) cos (ncot) dt = — < f(t) cos (ncot)\ l Jy /0 + nco / f (t) sin (ncot) dt > 

T { ’ J-T/2 J 

/((r/2)“)-/((-r/2) + )' cos (jicoT /2) + ncob n = ncob„ + (a + b) (—1)" = 0 


r T 11 
'-T/2 


on vérifie a’ n = ncob,, — ^cro cos (ncot 0 ) où cro est le saut à la discontinuité en to = k 

Ce qu’il faut retenir de cet exercice 

Vérifier que le comportement asymptotique des coefficients est compatible avec les propriétés 
de continuité de la fonction et de ses dérivées. C’est un bon test des calculs ; dans la question 
0 on a un exemple de série qui n’est pas dérivable terme à terme ; cela est à mettre en regard 
avec la non-convergence uniforme sur une période. 


1.6 On utilise 


© 


J e at sin (bt)dt 


ô—-r [-a sin (bt) 
+ b z 


bcos(bt)\ 


25 

















Chapitre 1 • Outils mathématiques de base 


et 


e cos (bt)dt = 


a 2 + b 2 

pour calculer les coefficients de Fourier. On obtient 

OO 


[b sin (bt) — a cos(èf)] 


1 O 

S(f ; t)= — (1 - e~ a ”) + — Y, [1 - (-l)"e 

7 TCI K ' TT ' L 


n= 1 


n- + a- 


COS lit 


La fonction / est continue donc S(f) = f Vf. En particulier 

/(O) = 1 = _(l- e —) + ^V [!-(-! Y‘e-^]- - 2 

tt a v tt z —' L J ?z z + a z 

«=1 

/(tt) = e -™ = — (i-«-“*) + —Y'I'i-c-i )"*-™ 1 (_1) " 

7 r « v ' tt L 


«=i 


n 2 + a 2 


Posons r,=f;-tLL et r, = 


« = 1 


72 2 + fl 2 


CO - 


n = l 


72 2 + fl 2 


1 = — (1 -e~ a U + — \Ti -e~ a7T TA 

ira x ' 77 L J 

= — (1 -e~ an ) + — \Ti -e~ a ”Ti] 

ira x ' 77 L J 


On en déduit 


T\ = 


7 T 


1 


2 a sinh( 7 ra) 2a 2 


T 2 = 


TT 


1 


2 a tanh(7ra) 2a 2 


D’après les notations de l’énoncé Si = — + T\ et S 2 = — + T 2 ; on a donc 

n L 


Sx = 


TT 


1 

+ ——x et S 2 — 


TT 


1 


2a sinh(7ra) 2a 2 “ 2 a tanh(7r«) 2a 2 

Ce qu’il faut retenir de cet exercice 

La valeur de la somme d’une série de Fourier en un point est donnée par le théorème de 
Dirichlet. Si la fonction est continue elle coïncide avec sa série de Fourier Vf. Les séries de 
Fourier fournissent un moyen efficace pour calculer la somme de séries numériques. 


1.7 Fa fonction est impaire donc a n = 0 et 

i r i r 

b n = — sin(û 2 0 f) sïn(nt)dt = — / [cos (72 — co Q ) t — cos (77 + co 0 ) t] dt 

TT . 2TT . 


2sin(n2ü7 t ) 72(-l)" ^ _ 2sinQo7r) ^ 


TT 


2 9 

(Oq — n z 


TT 


n =1 


2 7 

ù)q — n z 


sin (72 f) 
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Corrigés des exercices 


Si ct» 0 —1 alors s\r\(oj {] TT) -4 0 et tous les coefficients s’annulent sauf b\ qui vaut 

—2sm(ftj 0 7r) 

77 ("0 - !) _ 

Ce résultat est cohérent puisque la série de Fourier se réduit au seul terme sin t qui n’est autre que 
la fonction f{t). 

1.8 Le second membre de l’équation correspond à la fonction / de l’exercice 1.5 Q avec 
a = — 1 et ô = 1. On en déduit : 


b 1 = lim 
—> 1 


m 


tt 4 cos [(2k + 1) f] 

»+d 2 


Comme le second membre est une fonction paire et 2tt périodique et que l’équation préserve la 
parité (dérivées d’ordre 0 et 2), on cherche une solution sous la forme d’une fonction paire et 2tt 
périodique : y = ao/2 + a„ cos (nt). En injectant dans l’équation différentielle 

n 


y" — 4y = —4 (qq/2) + ( —n 2 — 4) a n cos (nt) 

n 

_ tt 4 cos [(2k + 1) t] 

2 “ ^ èS ( 2k + l) 2 

Par identification on obtient 

TT 4 1 

-T a 2k = 0 a ^k +1 = -TT-t-T 

4 - " TT (2k + if [(2k + if + 4] 

tt 4 ^4 cos [(2k + 1) t] 

~^ + ^h^ ik+i "> 2 [w+if+4] 

Remarquons que la série obtenue est bien deux fois dérivable puisque les coefficients de la série 

dérivée deux fois terme à terme se comportent encore en- - qui est le terme général d’une 

(2k d - 1 j 

série absolument convergente. 

Ce qu’il faut retenir de cet exercice 

La solution de l’équation est obtenue sous la forme d’une série rapidement convergente 
(a n ~ —) donc facilement calculable numériquement. 

FJ 4 


a 0 = 

y(t) = 


1.9 Pour avoir une série de sinus on construit un motif entre — L et L qui coïncide avec la 
© fonction / sur (0, L ) et avec —/ sur (—L, 0) ; on prolonge cette fonction impaire par périodicité 
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de période 2L. La série de Fourier de la fonction g ainsi obtenue ne comprendra donc que des 
termes en sin(/7 77—). Le coefficient de Fourier est donné par 


b n 


(2 L) J_ L 


/ g(x)sin (jitt— ^ dx = — f g(x)sin (jitt— ^ dx = — f f(x) sin ^mr —^ dx 

•L i/0 «/O 


Pour avoir une série de cosinus il faut construire un motif pair. 
• Série de Fourier cosinus de sin t sur [0, 7r| : 


2 f 77 4 1 

sin ? = flo/2 + cos 277? avec a n = — / sin ? cos 277? d? =-——- 

'' 77 /n TT An 1 — 1 

n> 1 


Série de Fourier sinus de cosh ? sur ]0,1 [ : 


/* 1 

cosh t = > h,, sin 7777 / avec b n = 2 / cosh? sin 7777? dt = Itt il — (— 1 )" cosh(l)l ——-- 

^ J n L J 77 2 77“ + 1 

n> 1 


Série de Fourier cosinus de cosh t sur ]0. I [ : 


cosh? = ao/2 + cosmrt avec a n = 2 / cosh ? cos 7? 77 ? d? = 2sinh(l)- 

n> 1 Jo 1 


(- 1 )” 
,2 77 2 + 1 


Ce qu’il faut retenir de cet exercice 

Une fonction peut être représentée sur un intervalle soit par une série de Fourier complète, soit 
par une série de sinus, soit par une série de cosinus. 


5 3 113 

1.10 On utilise la relation fonctionnelle : T(-) = --T(-) = - ypH, , 

Définition de la fonction B : B( 2, --) = F(2)F lIzlZS = r (2)r (l/2) / (-1/2) = _ 4 

’ 2 r(3/2) l/2r (1/2) 

De nouveau on utilise la relation fonctionnelle : 


277 ! 


277 ! 


2" (27?) (277 — 2) • • • (2) 2 2n n ! 


77 = 


77 (77 G N*) 


Pour calculer T(— 77 + -) on utilise la formule des compléments : 


r (77 h—) r(—77 h—) = —p-j—p - 

2 2 sin [77(77 + 2 )] 


1 2 2 ” 72 ' 

= (-l) n 77 =>- T(-n + -) = (-1)"— —'-yfrr (77 G N*) 

2 (277)! 
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Corrigés des exercices 


On effectue un changement de variable dans J n : 

1 


Jn 


poo 

/ f " É 
J o 


ln t dt = 


u=at 2 2\J~Ü J o \Cl 


u\ n / 2 du 

e (ln u — ln a) —= 
u 


1 


2 a 2 


r 


n + 1 \ , _, / « + 1 

— mn 1 


1.1 1 


1 jTT 1 

Jo = —[y + in(4 a)] J x = [y + ln(n)] 

2 \ a 2 a 

On fait le changement de variable u = fit. On obtient 

ricr + 1) 


h = 


2 <2+1 


On fait le changement de variable u = Vt + b 2 . On obtient 


fi = 


f e a(u h } 2 du = e ab Erfc (by/a) 
J b v=ufiâ V a 


© 


On identifie fi à la première représentation intégrale de la fonction B : 

2 

= 77 


1 1 




W 9 b 9 b“ 

ÿj On utilise l’identité : at~ + 2bt + c = a(t H —) + c - 

a a 


p OO -I p OO 

fi(x) = e ~ c+b2 / a / e- a(,+ ° )2 dt = e ~ c+bl/a -= / e^du 

Jx u=(r+%)fià \Ja 


2 e 


—c+b 2 /a 


Erfc 


+ X'Jci 


On identifie fi à la deuxième représentation intégrale de la fonction B : 


fi = 


l 


71 fi 1 1 + (X 1 — et 

(sin 0) a (cos 0)~ a dd = —B (—-—, —-—) = 


TT 12 


cos(7ra/2) 


0 Le cosinus intégral se comporte comme ln x pour x ~ 0 et il décroît exponentiellement 
à l’infini. L’intégrale existe donc. On intègre par parties : u = Ci(f) — > u' = co& ^ ; 


v' = cos {pt) —> v = — sin(pt) 
P 


fi = — sin(pt)Ci(f) 
P 


o P J o 


P J o 


cos (t) i r°° 

sm (pt)——dt = / 

t 2 p J 0 


sin [( p + l)f] + sin [( p — l)f ] 


dt 
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^ f°° sint TT f c 

On rappelle / - dt = — =>■ / 

Jo t 2 J g 


sin(û'f) , tt . 

—-- dt = — sign(a). Donc 


h 


TT 

4p 


[sign(p + l) + sign(p - 1)] 


+— pour p < — 1 

2p 

Opour \p\ < 1 
pour p > 1 

2p 

--pour p = ±1 


Cette fonction de p est représentée sur la figure 1.1. 



Figure 1.1 Intégrale 16 

^ On fait le changement de variable u = a (l + fi) ; il vient 

/■oo —u+ab 

h = - du = e ab E\(ab) 

Jab 11 

0 On fait le changement de variable u = a (t + b) ; il vient 

raiP+b) — u 

I % = e ab / - du 

Ja(a+b) 11 

Pour relier cette intégrale à une exponentielle intégrale il faut discuter du signe des bornes, 
sachant que la condition —fi ^ [a. b ] implique qu’elles sont toutes les deux de même signe 
(sinon u = 0 serait dans l’intervalle et l’intégrale serait divergente). 

Si 0 < a (a + fi) < a (fî + fi) — fi < a < f5 alors 


na(P+b) -U (• 

• ro o u 

f 00 e ~« ) 

h = e ab / - du = e ab \ 

/ - du — 

/ - d 11 f 

J a(a+b) M l 

. Ja(a+b) U 

J a(P+b) u J 


= e ab {E\ [a(a + fi)] — E\ [a(/3 + fi)]} 


Si a (a + fi) < a {fi + fi) < 0 a < p < — b on ne peut plus utiliser la relation de Chasles 
avec +oo car u = 0 serait dans l’intervalle ; on commence par faire un changement de variable 
pour se ramener à des bornes positives : 

pa(P+b) —u p-aiP+b) u 

h = e ab / - du = e ab / —du 

Ja(a+h) u u^-u J- a (a+b) 11 
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Corrigés des exercices 


Maintenant, à cause du e u , on utilise la relation de Chasles avec —oo pour assurer la conver¬ 
gence à l’infini. Ce faisant on introduit une divergence en « = 0 ; mais comme l’intégrale de 
départ est convergente on peut aussi la prendre en partie principale et on se ramène alors à la 
fonction E t : 


h 


ab 


r'—a(a+b) 


pp 


—du — PP 
u 


i'—a(f)+b) e u 

—du 

U 


e ab {Ei [a\a + b\] - Et [a\p + b\]} 


§) On fait le changement de variable u = ln t ; il vient 


h = 


ue 


—2m 


du f°° y/Ü 

— = / —z~e du 

U v=2u J q e~ u 


yfüe~ u du = T ( | j = 


l) On fait le changement de variable u = a(l + h) ; il vient 


ho — 


J ab 


sin(n — ab) 


du = 


J ab 


sin u cos ab — sin ab cos u 


du 


u 


cos (ab) ^ — — Si (ab)J + sin(ab)Ci(ab) 


On fait le changement de variable u = t 2n ; il vient 


/11 = 


ce» ,, 1 

U 2n 


-du 


Jab 1 + u 

On reconnaît la troisième représentation intégrale de la fonction B 

TT 


lu = — B(—, 1 - —) = — - - - 

2n 2n 2n 2n sin(7r/2/i) 


(n G N*) 


On fait le changement de variable u = t 2 /9 ; il vient 

In = l\^0-uŸT--Edu = i 2 
./() 2y/u 2 

81 5 1 243 

Y B( 2’2 “Tô" 77 


(1 - u) 3/2 u- 1/2 di. 


© 


Ce qu’il faut retenir de cet exercice 

Cet exercice met en évidence la diversité d’utilisation de ces fonctions définies par des inté¬ 
grales qui doivent désormais faire partie de la panoplie des fonctions usuelles au même titre 
que le sinus, l’exponentielle ou le logarithme. 
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1.12 D a) On a 

avec 


/ OO 

n 2 (t- 

-CO 

n ï(t-t') = 


t r )n 2 (t')dt l = J n 2 (t - 

1 pour t — \ <t' <t + \ 
0 sinon 


t')dt' 


Donc 


Intervalle en t 

Valeur de f* 2 (t ) 

f + l<—1 U f — 1 > 1 <F> f > 2 

-1 < t + 1 < 1 —2 < f < 0 

-1 < t - 1 < 1 0 < t < 2 

0 

J 

J 

/>/+1 

1 dt' = t + 2 

f dt' = 2-t 

t-1 


Finalement 

f* 2 (t) = (2 — \t\)u(2 — |f|) = 2A 2 (f) 

On peut utiliser une méthode graphique : n 2 (F) est un créneau fixe centré en F = 0 et 
ü 2 (f — t') un créneau « mobile » centré en t' = t. L’intégrale de convolution représente 
Faire du recouvrement des deux créneaux que l’on évalue en faisant varier t. 
b) Le calcul de l’intégrale ou la méthode graphique conduisent au résultat : 


r\t) 


n 2 (0 * 2A 2 (o 


= 0 

= ^ +3)2 
= 3 -t 2 

= ^- 3) 2 

= 0 


t < -3 
-3 < t < -1 
-1 < t < 1 
1 < t < 3 
t > 3 


c) Explicitons l’intégrale de convolution : 





t')g(t')dt' 




sin(t — t')(ci + t')dt' + 
— cos<7)sinf 



t')(a - 
- t')(a 


t'\)dt' 
t’)dt’ 


où le résultat s’obtient en intégrant par parties. 
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Corrigés des exercices 


d) Ces fonctions sont causales ; leur produit de convolution s’écrit donc 

fait)* f b (t)= . [ (t - t') a ~ l e~ (, ~ r ' ) t"’~ l e~ t ‘dt' 

r(fl)r(£) J o 


T(a)r(b) J o 

u(t)e~ f r ' 
T(a)T{b) 


f 

Jo 


(i t-t'f l t ,b x dt' = 




v =t'/t T(a)T(b) 


J 


(1 — v) a l v b l dv 


L’intégrale qui apparaît dans l’expression ci-dessus n’est autre qu’une fonction 

r(fz)r(z?) 

B(a,b) = — -—. En remplaçant on obtient : 


T(a + b) 


fa (0 * fb (t) = 


1 


Y(a + b) 


t a+b l e ' u(t) = f a+b (t) 


0 En utilisant la propriété de 8 {n) on a : 

r°° d n 

f (t) * s<»> (t ) = J fit - t')s (n V)dt = (- 1 )" — fit - 1 ') 

§) Produit de convolution de gaussiennes, 
a) On utilise l’identité :—at 2 + fit = —ait + fi/2a) 2 + fi 2 /4a 


= f M (t ) 


r'=0 


^dt = 


e -a( t+ f}/2a) e fJ /4a dt = g 


/3 2 /4a 


du 


' —OO 

n7 


u=y/cî(t+f$/2a) 


\ A 

a V 4 a 


b) Utilisons ce résultat pour évaluer l’intégrale de convolution de deux gaussiennes : 


fa it) * fb it ) = 


1 


2irab 


2iTctb 


,-'72 « 


e - (t - t 'fl2a 2 e -tV2b 2 dtl 

1 1 


exp 


—t' 2 |-h 

1 2a 2 2b 2 


exp 1 i?) *' 


_ /I 1 \ a 2 + b 2 

Cette intégrale est de la forme de celle de la question a) avec a = + 


et fi = t/n 2 . On a donc 

fa it) * fb it) = 


2 a 2 2b 2 ) 2 a 2 b 2 


1 


2TTClb 


,-r 2 /2a 2 


2ira 2 b 2 


» 2+i,! exp 4 « 4 èg 


2tt (a 2 + Z> 2 ) ^ 1 2 a 2 2a 2 (a 2 + Z? 2 ) 


b z t 


2 + 2 


2n t (a 2 + b 2 ) 


exp 


2 (a 2 + b 2 ) 


fcr(a,b )()0 


© 


OÙ (7 ici, b) = yj (fl 2 + Z> 2 ). 
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c) On a donc f* 2 (t) = f ay/5 ; à 1’ aide du résultat de la question précédente on calcule 

fa 3 = * fa = f aV 2 * fa = /^V) = faV 3 


On vérifie que ff n (t) = f a ^ (t) = ’- exp(--^). 

a v 277/7 2na- 

Ce qu’il faut retenir de cet exercice 

Les intégrales de convolution sont souvent difficiles à calculer ; on verra aux chapitres 2 et 3 
des méthodes plus efficaces que le calcul direct. Si les fonctions impliquées ont des formes 
simples - créneaux, triangles - la méthode graphique est efficace. 


Corrigésdesproblèmes 


Problème 1.1 

E X 

b„(t) sin(7i77y) ; en effet, cette série en sinus 


1 


satisfait automatiquement la condition aux bords. En supposant la série dérivable deux fois on 
obtient 


nfz 1 


£ w>- 


H7T 


— a + fî 


b n (t ) }> sin(«7Ty) = 0 


K(t) + 


n TT 


a + /3 b„(t) = 0 ^ b„(t) = B„ e~P'e-^f' 

u(x,t) = e~ pt ^ B„ e~ a( ^r> * sin(«77y) 
nfz 1 


Pour déterminer la suite de constantes arbitraires { B n } on utilise la condition initiale qui devient 
u(x, 0) = B n sin(«77y) = x(l — x) x G [0, /] 


n^l 


Les B n apparaissent donc comme étant les coefficients de Fourier sinus sur [0, /] de la fonction 
x{l — x). Donc 


B n = — / x(l — x) sin ( htt— ) dx = 


1 


4 1 1 1 - (- 1 )” 


7 T J 


Finalement 


4/ 2 v- 1 - (-1)" 

u(x,t) = — T e 13 Y -r- e ( ' ) si 

TT i n - 3 


n> 1 


sin(/î7Ty) 


On vérifie que la série est dérivable deux fois par rapport à x. 
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Corrigés des problèmes 


Ce qu’il faut retenir de ce problème 

La périodicité de la condition aux bords n(0, t) = u(l,t) = 0 suggère de chercher la 
solution sous la forme d’une série de Fourier. La valeur zéro de la fonction en x = 0 et 
x = / conduit à une série de sinus qui va automatiquement implémenter la condition. 


Problème 1.2 

On note /„ la fonction de période 1 qui coïncide avec B n sur ]0, 1[ et S„ = Sp(f n ) la série de 
Fourier de cette fonction. On notera {a k (n)\ et {b k (n )} les suites de coefficients de Fourier de 
cette série. 

f) On remarque que la fonction f\ est impaire, donc 


ffi(l) = 0 

r 1 î î 

b k ( 1) = 2 (x -)sin(27 Tkx)dx = - 

2 TTK 


La série s’exprime donc 


/o 


1 sin(27r^.x) 

1 7 T le 

k>\ 


Calculons le coefficient de Fourier a 0 (») de /« : 


pour n ^ 1 ao(n) = 2 f B n {x)dx = 2 f B' n+1 (x)dx = 2[B„+i(l) — B„ + i(0)] = 0 

Jq J o 

0 Calculons les autres coefficients de Fourier de f n : 

pour n ^ 2 : a*(n — 1) = 2 B n _ l (x)cos(2k7rx)dx = 2 / B' n (x)cos(2kTrx)dx 

J o J o 

= 2 B n (x)cos(2kTrx)\ l 0 + 2k7T x 2 / B n (x) sin(2kirx)dx 

J o 

= 2k7rbk(n) 

1 


On montre de même que 


bk(.n) = - — -a k (n - 1) 

Z77 k 


ak(n) = -——b k (n - 1) n > 2 
2jrk 


On en déduit que pour tout k ^ 1 


© 


a k (n) = - 


2tt k 


a k {n - 2) b k (n ) = - 


27 rk 


b k (n — 2) n ^ 3 
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On a des relations de récurrence par pas de deux entre les séries de coefficients. Ces relations 
sont initialisées par 


a k (1 ) = 0 

a k (2) = -^—b k { 1) = 


hW = - 


irk 


2irk 


2iT 2 k 2 


h( 2 ) = -^—a k (\) = 0 

2nTk 


Donc pour tout n impair la série des a k (n) est nulle ; de même pour tout n pair la série des b k (n) 
est nulle. Les autres coefficients s’obtiennent en résolvant les relations de récurrence, ce qui 
donne 


a k {2p) 


2(—1 ) p+1 1 2(-l^ +1 1 

(2ir) 7p k 2p k{ P (2 tt) 2 p + 1 k 2 P +l 


Finalement 

2(—l) p+1 ^ cos(2t Tkx) 2(-l )P +i ^ sin(27r*.r) 

Slp {X ~ (2 tt) 2 p 2^ k 2 P ~ p+l {X “ (2tt) 2 p +1 2s k 2 P + 1 

Je ^ 1 Je ^ 1 


Pour n ^ 2 les fonctions /„ sont continues du fait de la condition R„(l) = B n { 0). Donc 
f„(x ) = S n (x) = K n (x) Vx G [0, 1]. En particulier S n ( 0) = B n { 0). Pour n = 2 p on peut en 
déduire la valeur des fonctions de Riemann if (2 p) : 

S(2p) = E ^ = l(-V P+ \2ir) 2p B 2p (0) 

k^l * 


Grâce aux relations de récurrence entre les polynômes on trouve Bi{x) = -x 2 


2 X + 12 ’ 


1 


1 


1 


1 


TT~ 


77 


B$(x) = — x — —x + —x — On en déduit £(2) = —, £(4) = — 


24 


12 


24 


720 


6 


90 


Ce qu’il faut retenir de ce problème 

Les séries de Fourier fournissent ici un moyen de calculer la valeur de la fonction £ de 
Riemann pour des valeurs entières paires de la variable. 


Problème 1.3 

D O n va établir une représentation intégrale pour la somme partielle de Fourier. 
a) Exprimons la combinaison 


a n cos (nt) + b n sin(/7f) 


cos (nt)— / f (t 1 ) cos nt'dt' + sin(/7f)— / f(t')smnt'dt' 


77 


77 


— / f(t')cos(n(t' — t))dt' 

77 . 


36 











Dunod - La photocopie non autorisée est un délit 


Corrigés des problèmes 


© 


b) On en déduit que 

S N (t) = 


N i r fi N 

flo/2 + ^ a n cos (nt) + b n sin(«0 = — f(t') < - + ^ 

n^O 77 J -tt “ 

R>H* 5 


cos n(t' -() > dt' 


cos n{t' — t) > dt' 


On peut calculer explicitement la somme de cosinus en remarquant : 


N 

E 

n=0 


N 


cos nu = 


Rc ]T e in 


,i(N+l)u 


= Re 


- 1 _ Re { [ e i(N+ ^ u - l] [e~ iu - l] } 


n= 0 


e ,u - 1 


4 sin“(y) 


cos(Nu) — cos [(N + l)n] — cos (n) + 1 
4 sin 2 (|) 

2 sin(|) sin [(N + j)u] + 2 sin 2 (lj) 1 1 sin [(N + \_)u\ 


4 sin 2 (f ) 


2 2 sin (|) 


En reportant dans l’équation ci-dessus on en déduit l’expression pour la somme partielle 

S N (t) ■ 

1 r „.^sin [(N+ %)(?-t)\ 


S N (t) = — . 

2tt J_ n sin [4{t' - t)\ 


Appliquons cette représentation de la somme partielle à la fonction créneau 277-périodique 
impaire. 

a) Sa série de fourier n’a que des termes en sinus avec 


* = -!■ - 
77 


r 0 r 1 2 1 — (—i)" 

/ sin (nt)dt + / sin (nt)dt } = 


TT 


Onadonc/(0 = -V Sin(2 ^ + 1) ' 

TT ^ 


77 z —' 2k + \ 

k>0 


b) En utilisant le résultat de la question Q b) on a pour t G ]— 77 , 77 [ 


s»<o = A - 


sin [(« + !)((' -<)]., 


dt + 


sin [i N + \_)(t'-t)\ dt , 


v sin [\{t' - t )] J 0 sin [\{t' - t )] 

On pose u = t — t' dans la première intégrale et u = t' — t dans la seconde : 

sin(A + \)u 


S N (t ) = 


H-i; 


sin (N + \)u f n ' 
- r-=-—du + 


-du 


sin -u 


sin 2 u 
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Pour t G [0, 7 r] les intervalles [— t, tt — f] et [f, 77 + T] se recouvrent sur [f, 77 — f] ; sur cet 
intervalle la différence des deux intégrales s’annule ; il ne reste que 


1 sin UN + \)u\ 1 

S N {t) = — I - 1 2 J du - 


27 7 


sm 


(: l 2 u ) 


2tt 


sin [(N + \)u\ 
sin ( ju) 


du 


c) On regarde le voisinage de t = 0 ; en posant t = e les deux intégrales sont de la forme 


f 


g(t)dt = 2 eg(a) + o(e) 


sin UN + i)tl N 

où g(t) = --— -7-ê -—Dans la seconde intégrale a = tt et g( tt) = (—1) ; dans la 

sin( 5 t) 

première a = 0 et g(0) = 2N + 1. Donc dans la limite où N 
dominante. Si on pose m = N + -, v = mu on obtient 


oo la première intégrale est 


S N (t) ~- 


1 f m sin v 


77,77 _/ Q Sin (^) 


dv 


d) Pour obtenir le maximum : S' N (t) = 


1 sin,777 . , 

-qui s annule , 77 ? = 77 — > t M = 77 / 77 ? pour 

77sin(t/2j 

t < t M , sin(t/2) reste petit positif ; la dérivée est donc positive pour t < t M . On a bien un 
maximum de SV- 

v 

On calcule Sa?(?m = 77 /m) ; pour ,?7 ~ iV 1 on a — <S 1 puisque v G (0, 77 ) ; donc 

2m 


sm 


fJL) ~ JL 


V 2,77 


~ -— d’où 
2m 


çmax ^ 



On a .S// lx ~ —Si(77) ~ 1.18. Le fait remarquable est que cette valeur ne dépend pas de N. 

TT 

Quand N croît, ce premier maximum garde toujours la même valeur tandis que sa position 

îm — — se rapproche de zéro. 

N 

Ce qu’il faut retenir de ce problème 

Le phénomène de Gibbs traduit la convergence non uniforme de la série au voisinage d’une 
discontinuité. 
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Problème 1.4 

Q En utilisant les propriétés de la fonction T on a 


\nY(t)dt = / lnIXl — t')dt' = 


ln [r(f)r(l - t)] dt = - / ln 


2 
77 


Sin 77 1 - 


lnT(t)dt + / \nT{\ — t')dt' 


1 1 r 

dt = -\nir -/ ln(sin77t)^t 

2 2 ln 


- ln 77 + - ln 2 = ln VÏtt 
2 2 


Soit 


r x+\ 


a) On a pour, je > 0 


F(x) = / \nT{t)dt x ^ 0 


F'(x) = lnr(jc + 1) — lnr(.r) = ln — = hue 


n>0 


b) On en déduit pour jc > 0 


F(x) = xlnx — x + C 

»i 


On prolonge pour x —> 0 : F(0) = / lnr(f)c/t = ln \/277 = C =^ 

J o 

F{x) = xlnx — x + ln s/2tt pour x ^ 0 

0 Utilisons la formule des trapèzes pour F (x) 

r x+l i i c i 2 

F(x)= / lnT(t)dt = - [lnlX* + 1) + lnr(jc)]--lnr(jc) 

2 12 dx 


avec 8 G ]x, x + 1[ 


x=6 


Pour x grand 

On en déduit pour x grand 


1 d 2 

nJx 2 


ln r(.r) 


< 1 _1_ 

~ 12 jc — 1 


F(x) = — ln [T(x + 1 )r(x)] + o (—) = jc lnjc — jc + ln \/2tt 


On peut écrire 
1 


© 


2 ln [r(jc + l)r(jc)] = - ln 


r(jc +1) 


T(x + 1 ) 


= ln T(jc + 1) — ln \Jx 
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En remplaçant dans l’expression ci-dessus 

lnT(x + l) = x\nx — x + \n y/lïr + ln \[x + o(-) 

x 

on obtient que pour x grand (avec = 1 + o(l/x)) 

V(x + 1) ~ \f2Ïrx x+x ^ 2 e~ x [l + o(l/x)] 

Ce qu’il faut retenir de ce problème 

Cette démonstration de la formule de Stirling utilise largement les propriétés de la fonction 
Gamma. La formule de Stirling est plus connue comme approximation de la factorielle. 
En effet si on l’applique pour x = n G N : 

r(n + 1) = ni ~ V2n n n+1 e~ n 

Pour n = 5 on a 5! = 120 tandis que la formule de Stirling donne 118, soit moins que 2 % 
d’erreur. 
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